Elementary excitations and the phase transition in the bimodal Ising spin glass model 
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We show how the nature of the the phase transition in the two-dimensional bimodal Ising spin 
glass model can be understood in terms of elementary excitations. Although the energy gap with 
the ground state is expected to be 4J in the ferromagnetic phase, a gap 2J is in fact found if 
the finite lattice is wound around a cylinder of odd circumference L. This 2J gap is really a 
finite size effect that should not occur in the thermodynamic limit of the ferromagnet. The spatial 
influence of the frustration must be limited and not wrap around the system if L is large enough. In 
essence, the absence of 2 J excitations defines the ferromagnetic phase without recourse to calculating 
magnetisation or investigating the system response to domain wall defects. This study directly 
investigates the response to temperature. We also estimate the defect concentration where the 
phase transition to the spin glass glass state occurs. The value pc = 0.1045(11) is in reasonable 
agreement with the literature. 

PACS numbers: 64.60. ah, 75.10.Hk, 75.10.Nr, 75.40.Mg 



I. INTRODUCTION 

Spin glassesi"— have attracted much interest for quite a 
while. Due to the considerable complexity of real materi- 
als much computational effort has gone into studies of a 
simplified model'"' that is nevertheless thought to include 
the essential ingredients that lead to spin glass behaviour. 
Since even this model is not trivial, a considerable indus- 
try has developed over time devoted to particular models 
that, although probably unphysical, have provided sub- 
jects for the development of numerical techniques^'*'. 

Systems known as spin glasses are disordered magnetic 
systems characterised by a random mixture of ferromag- 
netic and antiferromagnetic exchange interactions lead- 
ing to frustration^. Typically, at low temperatures below 
a critical temperature Tc, a system undergoes a phase 
transition from a ferromagnet to a spin glass at some crit- 
ical concentration pc of antiferromagnetic interactions. 

The model studied in this work is the bimodal, or 
± J, Ising spin glass in two dimensions. This system has 
quenched bond (short range, nearest neighbour) interac- 
tions of fixed magnitude J but random sign. The con- 
centration p of negative, or antiferromagnetic, bonds is 
varied from zero up the canonical spin glass aX p — 0.5. 
It is believed that the spin glass can only exist at zero 
temperature^ where p > pc withi^iii pc = 0.103. This 
is clearly below the concentration p„ = 0.109 at the (fi- 
nite temperature) Nishimori point, indicating a reentrant 
phase transition as confirmed by Monte Carlo workiS. 

The ground state is highly degenerate with an entropy 
per spinHiii of 0.07fc. Consequently spin correlation 
functions are not guaranteed to take values ±1.0. If a 
nearest-neighbour bond correlation function does have a 
value ±1.0 then we call that bond a rigid bond^^. This 
means that the spin alignment across the bond is the 
same in all ground state configurations. A recent studyiS 



suggests that the rigid lattice does not percolate in the 
spin glass phase. This is consistent with the idea that 
the ferromagnetic phase is characterised by percolation 
of rigid bonds. 

Droplet theoryi^r— has enjoyed much success with re- 
gard to understanding the spin glass phase. The essential 
idea is that reversing all the spins in a compact cluster 
with respect to a ground state provides a low energy exci- 
tation. Typical droplet excitations dominate the thermo- 
dynamic behaviour. A closely related idea is the domain 
wall defect^^'^'^; essentially a droplet perimeter that ex- 
tends to infinity. With a continuous distribution of disor- 
der these related views seem to be equivalent^- according 
to the predictions of droplet theory. 

For the bimodal model domain wall defects have, for 
example, been appliediii^ to the determination of the 
value of the critical defect concentration pc- Neverthe- 
less, it still remains unclear whether droplet theory is 
appropriate^^. The ground state is not unique and a 
droplet may represent some different ground state; not 
an excitation. 

For this study the Lx L square lattice is wound around 
a cylinder, that is we use periodic boundary conditions 
in one dimension. In the second dimension the system 
is nested in an infinite unfrustrated environment. There 
are no open boundaries. If the circumference L of the 
cylinder is even then the energy gap is 4 J. Otherwise it 
is 2 J. In the spin glass phase the distribution of degen- 
eracies of the first excited state is extreme with a long 
tail representing large values We have also looked at 
systems with open boundaries and have found extreme 
distributions of 2 J excitations in some agreement with 
Wang^^. 

The issue of the size of the energy gap of the bimodal 
Ising spin glass dates back to the proposal of Wang and 
Swendser>^ that it should be 2 J in the thermodynamic 
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limit. It now seems clear that in fact there is no energy 
gap at all and the low temperature specific heat varies as 
a power law ^ T^". The first indications of this ap- 
peared in Ref. 31 and confirmation'^^ from the evaluation 
of very large Pfaffians has recently appeared. 

The issue that remains unclear is the value of the crit- 
ical exponent a. For the case of continuous (Gaussian) 
disorder, direct calculations'^-'-- report that the specific 
heat is linear with a — — I.O. For bimodal disorder 
Monte Carlo work^ reports that a = —4.21 while droplet 
theory2E suggests that a = —3.0 although the tempera- 
ture range used is extremely narrow. Other Monte Carlo 
results^i^ for the correlation length with the assump- 
tion of hyperscaling gives a = —7.1. Universality is hard 
to prove. 

The exponent a is difficult to estimate. One reason 
that makes this so for the bimodal case is that the spe- 
cific heat is not normally distributed. We have performed 
some calculations with open boundaries and find that the 
distribution of the specific heat has a tail for low tem- 
perature and small values of linear sample size L. The 
methods used were direct evaluation of Pfaffians as well 
as summing the density of states^. Although it is rea- 
sonable to believe that the specific heat will be normally 
distributed in the thermodynamic limit, it is not clear 
what value to use from calculations with finite L. 

It is at least clear now that the low-temperature spe- 
cific heat contains contributions from excitations having 
a range of energies. This fits well with droplet theory'^^ 
where it is predicted that a = 1 — 2/9s with the frac- 
tal dimension of domain walls given by = 29s. If 
0s = 0.5 as reported^i^^i^ then a = -3.0. How- 
ever, other work"^^"^^ predicts values d/ > 1 that imply 
a > —3.0. It seems unlikely that droplet theory can pre- 
dict a value in agreement with a = —4.21 or a = —7.1. 

To obtain a simple description of the ferromagnetic 
phase we can start with the case of low defect concen- 
tration p. The defect bonds are widely separated and 
the ground state is unique (aside from global inversion). 
So the degeneracy of the ground state is Mq = 1. We 
can find first excited states by flipping a spin at either 
end of a defect bond. Thus the degeneracy of the first 
excited state is Mi — ApN where the square lattice has 
N sites and 2N bonds. The value of the density of states 
per spin is 4p. We can think of clusters of disorder 
each composed of one negative bond and two frustrated 
plaquettes. 

As the concentration p < pc increases the clusters of 
disorder grow in size and influence. The distribution of 
the density of states becomes less normal and its peak 
moves above 4p. Nevertheless, the rigid lattice still perco- 
lates and there remains some finite magnetisation. With 
a lattice of finite size, wound in one direction, the 4J 
excitations occur in two classes. Some are derived lo- 
cally and are not influenced by the boundary condition; 
just like the simple case of low concentration. Others are 
formed by extending all the way around the system. 

Since it is not easy to distinguish between these two 



classes, we employ the device of flxing an odd value of 
the circumference L of the cylinder. In this case the 2 J 
excitations are entirely nonlocal. Fig. [1] shows an exam- 
ple. The excitation depends on the boundary condition 
and would not exist otherwise. All 2 J excitations involve 
flipping all spins on one side of some closed path around 
the system. Other closed paths can give excitations with 
energies equal to an odd multiple of 2 J. A 4 J excitation 
requires two paths. 




FIG. 1: (Color online) An example of a 2J excitation. On 
the left is a ground state configuration with six frustrated 
plaquettes and five unsatisfied (jagged) bonds. On the right 
is a first excited state obtained by flipping all spins on one 
side of the vertical broken fine. The excited state has six 
unsatisfied bonds. Periodic BCs are indicated by the top and 
bottom dashed vertical lines. 

Our main message here is that it is possible to essen- 
tially define the ferromagnetic phase by the absence of 
these 2J excitations. Alternative approaches'^ include 
the imposition of domain wall defects and the calcula- 
tion of magnetisation. These nevertheless lack clear sys- 
tematics due to the large degeneracy of the ground state. 
Domain wall defects may not represent excitations at all 
since they can correspond to alternative ground states. 
Sampling of domain walls cannot be done in a controlled 
way and it is not obvious^^ how we can obtain typical 
representative domain walls. 

Calculation of the magnetisation is also problematic as 
a result of the ground state degeneracy. In Ref. m for 
example, the algorithm starts with a ground state and 
proceeds with a Monte Carlo simulation to determine a 
typical value of the magnetisation. 

In this work we propose a simple picture of the ferro- 
magnetic phase that is evaluated from the response to 
temperature alone. The number of lowest energy exci- 
tations is counted exactly. In the thermodynamic limit 
these excitations can only exist in the spin glass phase. 
Details of our results are given in Sec. IIIII after a brief 
account of our method. 



II. FORMALISM 

We use the Pfaffian method and degenerate state per- 
turbation theory to calculate the degeneracies of the ex- 
cited states. The planar Ising model can be mapped 
onto a system of noninteracting fermions. Each bond 
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is decorated with two fermions, one either side. A square 
plaquette then has four fermions inside and four others 
across the bonds, as shown in Fig. [2] For a system with 
N lattice sites we have AN fermions in total. The parti- 
tion function can be written a.a^^d^ 
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At arbitrary low temperature the internal energy can 
be expanded as^^ 



m=l 



(5) 



where the coefficient Um is expressed as 

Urn = -2" J Tr i?" (6) 

with 



o 



FIG. 2: A plaquette with associated lattice fermions. The 
filled circles are associated with the plaquette, and the pairs 
of filled and open circles are associated with the bonds. 



Jl cosh( Jy/fcT) I {detD^/^. (I) 



The product is over all nearest-neighbor bonds Jij on 
an A'' site lattice. The matrix £) is a 4A x AN 
skew-symmetric matrix that comprises constant diago- 
nal blocks, and off-diagonal blocks that depend on tem- 
perature T through matrix elements ± tanh Jij/fcT. The 
factor (detZ))^/^ is precisely the Pfafhan^i^. This for- 
malism is applicable to any distribution of disorder. For 



the bimodal model J,, 



±J. 



At zero temperature there are defect eigenstates of D 
with eigenvalues equal to zero. Each defect eigenstate 
can be expressed as a linear combination of the fermions 
localized in a frustrated plaquette. The number of these 
defect eigenstates is exactly equal to the number of frus- 
trated plaquettes. At low temperature each defect eigen- 
value approches zero as 



e — ±—X exp 

2 ^\ kT 



(2) 



where r is an integer and A" is a real number. These 
quantities r and X can be obtained using degenerate 
state perturbation theory^. The ground state energy 
is written as 



Uo = -2NJ + 2jJ2rd 



(3) 



where the sum are over all defect eigenstate pairs. The 
ground state degeneracy is 



(4) 



and the ground state entropy can then be written as 5*0 



R = Dig,i{l + DiGi){l + D2G2) • • ■ (1 + i?r„_a„„J. 

(7) 

The 2x2 block diagonal matrix Di is defined accord- 
ing to D = Dq + SDi where Dq is the matrix D when 
T = and 6=1 — tanh J/kT. Di has non-zero ma- 
trix elements joining two fermions across bonds only. 
The 4x4 block diagonal matrix Qci is derived from 
the continuum Green's functior*^ and has matrix ele- 
ments connecting the fermions within a plaquette. D2 
is given by D2 = DigdDi and, for r > 2, Dr = 
Dr-iil + Gr-2Dr-2) ' ' ' (1 + GiDi)gciDi. The Green's 
function Gr is given by*^ 



N(r) 

i=l ^ 



(8) 



where \r,i) is the ground state defect eigenstate of D 
with eigenvalue e* . The integer r represents the order of 
perturbation theory at which the degeneracy is lifted. It 
is also the index r in Eq. The total number of such 
eigenstates is N{r). 

The coefficient Um can also be expressed in terms of 
the degeneracies of the excited states. We denote the 
degeneracy of the zth excited state as Mi. The partition 
function of the bimodal Ising model can be expressed in 
terms of the degeneracies as 



"^0 

Z = 2Moe^-pT ( 1 



Ml _2J_ M2 _4J 

e ''^ H e >'T 

Mo Mo 



(9) 



Using some thermodynamic relations together with the 
expansion of InZ using the Taylor series ln(l -\- x) = 



., we obtain, for example. 



U2 
U. 




(10) 



From these relations the ratios -jf- for all excited states 
can be obtained recursively. Note that M\ is the number 
of 2 J excitations. 
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III. RESULTS 

We have calculated ^ for system sizes up to L = 129 
and concentrations p ranging from 0.050 to 0.150. The 
number of disorder realizations ranges from 20000 for the 
smallest size to 2000 for the largest. We denote as Pi 
the probability of finding ^ > 0. In Fig. H Pi is 
plotted as a function of system size L for various defect 
concentrations. The error bars are evaluated using the 
bootstrap method^-S. The transition concentration pc is 
indicated where the L dependency of Pi changes from de- 
creasing to increasing. We can see that Pi is decreasing 
for p < 0.102. The system can be regarded as ferromag- 
netic below this concentration. Since Pi is increasing 
with L ioT p > 0.106, the system can be regarded as a 
spin glass. We conclude from these results that the value 
of Pc lies between 0.102 and 0.106. 




FIG. 3: (Color online) The probability Pi of finding > 0, 
plotted as a function of system size L for various values of 
antiferromagnetic bond concentration p. 

We have done a scaling plot using the relatior^, 

PiL^ ^ fiip-Pc)L'f'). (11) 

It is reasonable to fix ^ = since the value of Pi is 
bounded to the range [0,1]. In any case, with ^ not fixed, 
the best scaling plots have tp < 0.001. The parameters pc 
and (j) are chosen to minimize the quality paramete r^^i^'' 
S. The best fits give pc = 0.1045(11) and ^ = 0.532(72) 
with S = 0.62. The resulting scaling plot is shown in 
Fig. m The error bars of each parameter are obtained 
using the method described in Ref. [Til We fix the cor- 
responding parameter at various values and minimize S 
with respect to the other parameter. The range of the 
fixed parameter that gives S double the minimum value 
is regarded as the error bar. For example we show in Fig. 
[5] the variation of the partial minimized value of S(j)c, <p) 
as a function of Pc- The error bar of (p can be obtained 
in the same way. 

The above value of Pc agrees, within error bars, with 
Pc = 0.103(1) proposed in Ref. [111. Note that we also 



performed the analysis using data from systems with L < 
65 and get Pc ^ 0.105. This indicates that the effect of 
finite size is the overestimation of Pc- It is expected that 
if we perform this analysis using data with L > 129, we 
will get a smaller value of Pc- 
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FIG. 4: (Color online) The scaling plot of Pi as a func- 
tion of the antiferromagnetic bond concentration p with pc = 
0.1045(11) and (j) = 0.532(72). 
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FIG. 5: The variation of Smin as a function of pc. 

We have also investigated the distributions of the 2 J 
excitations in the ferromagnetic phase. We denote as 
Ci{x) the probability of finding < x. In Fig. [HI 
Ci{x) with p — 0.090 is plotted for various values of 
L. It is clear that the most likely value of is zero. 
The probability of getting ^ > is decreasing with L. 
We may expect that in the ferromagnetic phase the 2 J 
excitations vanish in the thermodynamic limit. 

Since there are no 2 J excitations in the ferromagnetic 
phase in the thermodynamic limit, the first excited state 
has energy 4J. We have investigated the behavior of 
the 4 J excitations by calculating the ratio ^ for system 
sizes up to L = 97. We denote as H2{x) the probability 
density function of getting jji^^ = x. We use the ker- 
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FIG. 6: (Color online) The probability Ci{x) of finding ^ < 
a; for p = 0.090. 
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FIG. 8: (Color online) The probability density H2{x) of get- 
1 Ma 
Mo 



ting -X^ = X ioT p = 0.090 with even L. 



nel density estimation algorithm^ to obtain H2{x). In 
Fig. [71 H2{x) with p = 0.090 is plotted for various odd 
values of L. A sharp peak develops with increasing L. 
We expect to get a definite value of in the ther- 

modynamic limit. It is interesting that tliis behaviour 
does not depend on whether L is odd or even. In Fig. HI 
H2{x) with p — 0.090 is plotted for various even values 
of L. The distributions are much the same and provide 
the same conclusions. From these results we have that 
the energy gap in the ferromagnetic phase is 4 J. 

We can expect this also from the behavior of the spe- 
cific heat at low temperature. When the temperature is 
low enough the behavior is dominated by the first excited 
state and can be expressed as^ 
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FIG. 7: (Color online) The probability density H2{x) of get- 
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Sharpness of the distribution of t^-t^ satisfies the re- 



that in the ferromagnetic phase the specific heat has a 
definite value in the thermodynamic limit. At low tem- 
perature Cy is proportional to exp(— 4J/fcr) and the en- 
ergy gap can be regarded as 4 J. 

The distribution of the 2 J excitations in the spin glass 
phase is quite different. In Fig. [HI Ci{x) withp = 0.110 is 
plotted for various values of L. Although the most likely 
value of ^ is still at zero, the probability of getting 

4t" > is increasing with L. The distributions of 41^ do 

Ado ° Ma 

not have a sharp peak but broaden when L is increasing. 
We have that the 2 J excitations persist as L increases. 
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FIG. 9: (Color online) The probability C\{x) of finding 
X for p = 0.110. 
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quirement of as a physical quantity. We can expect 



The distribution of the 4 J excitations in the spin glass 
phase is also different from that in the ferromagnetic 
phase. In Fig. [101 H2{x) with p = 0.110 is plotted for 
various odd values of L. The most likely value of H2{x) 
increases with L and the distributions broaden. This be- 
havior of H2 {x) is similar— to that of the canonical spin 
glass {p = 0.5) with even L. In particular, the height 
of the peak of the distribution collapses with increasing 
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L. We have also checked the distributions of H2{x) with 
p = 0.110 and even L. The results are shown in Fig. [TT] 
The characteristics are the same for both odd and even 
L. 
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FIG. 10: (Color online) The probability density H2{x) of get- 
ting = j: for p = 0.110 with odd L. 
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FIG. 11: (Color online) The probability density H2(x) of get- 
ting T^-fj^- = X for p — 0.110 with even L. 



romagnetic phase is characterised by the absence of low- 
est energy, that is 2J excitations. Our method counts 
the number of excitations exactly without bias. It is not 
necessary to work with some typical ground or excited 
state. 

Distributions of the number of 2 J excitations are 
shown to differ in character between the phases. In the 
ferromagnetic phase the number declines as the (odd) 
circumference L of the cylindrical winding increases. A 
finite-size scaling analysis produces a data collapse of ex- 
cellent quality to support our conclusion that 2 J exci- 
tations do not exist in the thermodynamic limit of the 
ferromagnetic phase. In the spin glass phase the situa- 
tion is reversed with the degeneracy of the first excited 
state increasing with L. 

The energy gap in the ferromagnetic phase is 4 J. For 
even values of L the first excitations have energy 4J. 
We have also presented distributions of 4J excitations 
so as to indicate that there is no essential dependence on 
whether L is even or odd. In the ferromagnetic phase the 
peak grows taller and narrower with increasing L and will 
presumably lead to a unique value of the low temperature 
specific heat in the thermodynamic limit. 

In the spin glass phase the behaviour of the distribu- 
tions is quite different. Essentially they are extreme with 
long tails. As L increases the tails become fatter and the 
peak collapses. We believe that this is consistent with 
a power law behaviour for the low temperature specific 
heat. The extreme distributions only indicate spin glass 
behaviour; a proper statistical mechanical description of 
the model requires a summation over the entire density 
of states. This seems to suggest that thermally active 
droplets can indeed take many different values of energy. 

Finally, we have not found any evidence that indicates 
a random antiphase statc"'^^, although we cannot rule out 
a situation where percolation of rigid bonds coexists with 
zero magnetization. 
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